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Preface 


The purpose of this study was to model the effect charged particle radiation 
has on the working parameters of metal-oxide-semiconductor field-effect transis¬ 
tors (MOSFETs). There are two individuals I would especially like to thank. One 
is my thesis advisor Michael Sabochick, for the helpful advice and guidance he 
provided throughout the project. The other is Roger Tallon of the Air Force 
Weapons Laboratory, who originally suggested this as a possible thesis topic, and 
provided an (at the time) unpublished copy of their experimental results. 
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Abstract 

A simple model to characterize the effects of heavy charged particles on 
metal-oxide semiconductor field-effect transistors (MOSFET) is presented. The 
model is applied to experimental results provided by the Air Force Weapons 
Laboratory, and an attempt made to simulate saturation phenomena observed in 
the threshold voltage change. The model assumes all trapped holes are within a 
few nanometers of the oxide-semiconductor interface, and takes into account the 
resultant counter electric field, and its effect on the yield fraction escaping recom¬ 
bination. An equation relating threshold voltage change as a function of dose is 
derived and used, but does not duplicate the saturation effect. This is because 
charge trapped near the interface reduces the internal oxide electric field very lit¬ 
tle. However, similarities in p-channel and n-channel MOSFET damage sensitivi¬ 
ties indicate considerable bulk charge trapping. This would modulate the internal 
electric field considerably. To determine whether this accounts for saturation will 
necessitate refinement of the model, taking hole transport and bulk trapping into 


account. 








MODELING THE EFFECTS OF HEAVY 


CHARGED PARTICLES ON METAL-OXIDE- 
SEMICONDUCTOR FIELD EFFECT TRANSISTORS 

I. Introduction and Background 

One of the most important devices used in very-large-scale integrated circuits 
is the metal-oxide-semiconductor field-effect transistor (MOSFET). Figure 1 
shows a schematic of an n-channel MOSFET. The device works as follows. 

When a positive voltage is applied to the gate, a strong electric field sets up 
through the gate oxide and penetrates into the p-type Si substrate. In the region 
just beneath the gate oxide, this field repels the majority carriers (which, for n- 
channel MOSFETs, are holes) and attracts minority carriers (electrons), thus 
forming an inversion layer. With a potential difference applied between source 
and drain, current flows through this inversion layer and the device is considered 
"on". The operation of a p-channel MOSFET is similar, only the p-type and n- 
type materials are reversed, and a negative voltage is applied to the gate contact. 

MOSFETs are very sensitive to the effects of ionizing radiation. The 
sequence of events associated with irradiation of an n-channel MOSFET is shown 
in Figure 2. Ionizing radiation, such as gamma-rays, electrons, and heavy charged 
particles, initially create many electron-hole pairs throughout the gate oxide. A 
certain fraction of these pairs are then annihilated through recombination. The 
positive applied voltage rapidly sweeps the remaining electrons from the oxide and 
into the metal gate contact. The holes transport more slowly toward the gate- 
semiconductor interface, where a certain fraction are caught in deep hole traps 
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Figure 1. Schematic of an N-channel MOSFET (Ref 1:22) 

and the rest disappear into the semiconductor. These trapped holes permanently 
alter the electric parameters of the device, primarily by lowering the threshold 
gate voltage required to create the inversion layer. In fact, if enough holes are 
trapped at the interface, an inversion layer may form without applying any vol¬ 
tage at all. For p-channel MOSFETs, the process is similar, except the holes 
migrate toward the metal-gate contact and are trapped at that interface. 

Recent experiments have been conducted by the Air Force Weapons Labora¬ 
tory to characterize the effects of ionizing radiation on MOSFETs. References 3 
and 4 describe this work in detail. Large numbers of non-radiation-hardened 
Intersil 3N161 p-channel and 3N171 n-channel discrete MOSFETs were irradiated 
with protons, electrons, and Co-60 gamma rays. .After exposure to specified levels 
of radiation (under various bias conditions), the transistor gate threshold voltages 
were measured "in situ", and the changes in threshold voltage found as a function 
of radiation angle of incidence (the angle between the proton track and the electric 
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Figure 2. Sequence of Events Associated with Irradiation of MOSFET (Ref 
2:67) 


field), applied field, particle energy, and total ionizing dose. 

Of particular interest were the results of the proton irradiations. In this case, 
the MOSFETs were irradiated at room temperature with protons of energies from 
2-16 MeV. The angle between the proton track and the electric field was varied, 
and data collected for the cases of 0, 45, and 80 degrees. Figures 3 and 4 show 
the damage sensitivity (change in threshold voltage per dose) for the p-channel 
and n-channel MOSFETs, as a function of particle energy and angle of incidence. 
In general, the damage increases with both particle energy and angle of incidence. 
This is as expected, as will be explained in the theory section of this paper. 

What is not expected, however, is the anomaly where the damage caused by 
the 45 degree protons is greater than that of the 80 degree protons for particle 
energies less than 6 MeV. When first observed, the researchers thought that this 
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anomaly was due to the additional material (silicon dioxide passivation and 
aluminum metallization) the 80 degree protons had to pass through before reach¬ 
ing the gate oxide. If these additional material lengths were greater than the 
mean proton ranges, then the number of 80 degree protons reaching the gate oxide 
would be reduced. However, the actual distance traveled by the protons was 
found to be 6 to 7 times less than the mean ranges of 2-6 MeV protons. This 
implied that all of the protons reached and passed through the gate oxide, regard¬ 
less of the angle of incidence. 

In a further attempt to understand this anomaly, the data was replotted in a 
different format, shown in Figures 5, 6, and 7. Figure 5 shows the change in 
threshold voltage as a function of dose for the various angles of incidence. For 
this case, the gate bias potential was -5 Volts. Figures 6 and 7 show essentially 




























Figure 4. N-channel MOSFET Damage Sensitivity as a Function of Particle 
Energy and Angle of Incidence (Ref 4:1209) 


the same thing, though with different gate bias potentials. On the basis of these 
graphs, the researchers deduced that a saturation effect was occuring within the 
test sample, that had the same effect as reducing the internal electric field during 


To explain this phenomena, the researchers proposed that the buildup of 
trapped holes at the interface resulted in a counter electric field. The yield frac¬ 
tion of holes escaping initial recombination increases with both the angle of 
incidence and the magnitude of the electric field. Initially, the higher yield frac¬ 
tion of the 80 degree protons results in a larger charge buildup at the interface. 
However, this larger buildup offsets the electric field within the oxide, decreasing 
the fraction of holes escaping recombination. Thus, the threshold voltage change 
due to 80 degree protons saturates before the 45 degree case. 
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Figure 5. Change in Threshold Voltage as a Function of Dose and Angle of In¬ 
cidence for Applied Gate Voltage of -5 Volts (Ref 4:1212) 


The goal of this thesis project was to give this rough qualitative model a 
more precise quantitative formulation, and verify that the counterfield effect actu¬ 
ally explains the anomaly. The next section describes the stages shown in Figure 
2 in more detail. Assumptions used in the model are stated, and mathematical 
formula presented. The following section presents the major results and analysis, 
and is, in turn, followed by the conclusion. 
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Change in Threshold Voltage as a Function of Dose for 80 Degree 
Angle of Incidence and Gate Voltages of -5 and -10 Volts (Ref -4:1213) 
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Change in Threshold Voltage as a Function of Dose for 45 Degree 
Angle of Incidence and Gate Voltages of -2 and -5 Volts (Ref 4:1213) 







II. Theory 


This section describes the assumptions and mathematical equations that 
enter into each stage of the model. Considered in turn are columnar phenomena 
including initial electron-hole formation and recombination, bulk phenomena 
including electron and hole transport, and interface phenomena including deep- 
hole trapping, and radiation-induced traps. Following this, an equation relating 
the change in threshold voltage to the total ionizing dose is developed. Both the 
dependence of the electric field on trapped charge, and the dependence of yield 
fraction escaping recombination on electric field are considered in this develop¬ 
ment. 


Columnar Phenomena 

As the protons pass through the oxide, they lose energy through ionization 
and excitation of electrons across the energy bandgap; electron-hole pairs are 
formed. The amount of energy needed to create an electron-hole pair in silicon 
dioxide has been estimated to be about 18 eV (Ref 5:1520). Protons, as well as 
other heavy charged particles, leave a very dense column of electron-hole pairs in 
their track. However, to determine the initial charge produced by an increments 
dose, the electron-hole pairs are assumed to be evenly distributed throughout the 
bulk of the oxide. This incremental charge is given as 

AC = K g ADAlf v (1) 

where 


AC = charge (Coulomb) 





Kg = generation constant (Coulomb/cm 3 —rad) 

AD = dose (rad—Si0 2 ) 

A = area of oxide (cm 2 ) 

1 = length of oxide (cm) 
f y = fraction yield 

The generation constant, Kg, is found from 

„ _ (.OlJ/kg—rad)(1.609 X 10~ 19 C/hole)p 
** (1.609 X 10" 19 )(1000 g/kg)W 


p = density (g/cm 3 ) 

W = energy per electron—hole pair formed (eV/hole) 

For Si0 2 the density is 2.2 g/cm 3 and the energy required to form an electron hole 
pair is 18 eV. Filling these numbers into Eq (2) gives a generation constant of 
1.222 X 10 -6 Coulomb/cm 3 —rad. 

Immediately after the electron-hole pairs are produced in the dense column, a 
certain fraction of them will recombine. The fraction that survive. f v , is a func¬ 
tion of the particle type and energy, the electric field within the oxide, and the 
incident angle between the particle track and electric field. To find the fraction 
yield, a numerical procedure developed by Oldham (described in references 6 and 
7) is used. A description of this procedure follows. 

Figure 8 shows a schematic of the track of a charged particle after passing 
through the silicon dioxide insulator. Looking down the direction of the track, 
the initial charge density configuration is assumed to be 
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Figure 8. Charged Particle Track Through the Oxide 


n ± = electron—hole density (cm -3 ) 

N 0 = initial linear density (cm -1 ) 
b = Gaussian radius (cm) 
r = Distance center of the track (cm) 

The electrons and holes then move according to the equation (Ref 7:2695) 

dn. dn. 

= D ± V 2 n ± + n ± E x -^ -an + n_ 


D ± = Diffusion coefficient (cm 2, s) 

H ± = mobility (cm^'V-s) 

E x = Component of electric field in the x—direction (V/cm) 
a = Recombination coefficient (cm 3/ s) 

The first term on the right hand of Eq (d) represents diffusion, the second term 
represents drift under the influence of the external field, and the third term 
represents recombination. 
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The actual program (found in the Appendix) sets up a grid in the x and in¬ 
direction with z in the direction of the proton track, and y perpendicular to the 
electric field for all proton angles of incidence. The initial electron-hole density is 
found from Eq (3) for each point in the grid. These densities are then placed into 
a finite difference form of the right hand side of Eq (4), and dn ± /dt found at each 
point on the grid. The new densities are then found from 


where At is the increment in time. These new densities are inserted into the right 
hand side of Eq (4), and the process repeated until the electron and hole columns 
separate under the influence of the external field. Integrating n + (r,t) over the area 
of the grid gives a new linear hole density N(t). 

Figure 9 shows the electron-hole columns moving past each other under the 
influence of an electric field of arbitrary direction. For the case of a particle track 
perpendicular to the electric field, the fraction yield, f y , will be given simply by 
N(t)/N 0 , where N(t) is the linear density after the columns have separated. The 
analysis up to this point is only good for the region between the dashed lines in 
Figure 9, but it can be generalized to three dimensions. If it is assumed that the 
hole column to the right of the right-hand dashed line remains unchanged, then 
the fraction yield can be found from 


where T is the time required for the right hand of the electron column to cross the 
left hand of the hole column. 

The parameters n ± , D ± , a, N 0 , and b in Eqs (3) and (4) are found in the fol¬ 
lowing way. The room temperature mobilities of silicon dioxide are found quite 
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Figure 9. Electron and Hole Columns Moving Under the Influence of Arbitrary 
Electric Field (Ref 7:2696) 


readily from the literature. The values used in this analysis were 
H+ - 2 X 10 -5 cm 2/ V—sec (Ref 8:2012) and - 20 cm 2/ V-sec (Ref 9:1333). The 
diffusion coefficients, D ± , are then found from the mobilities by using the Einstein 
relation (Ref 7:2696) 




M ± kT 


where k is the Boltzmann constant. T is the temperature, and e is the electron 
charge. The recombination coefficient, a. is found by using an expression derived 
by Langevin (Ref 7:2696) 

(/i+ + H.)e 

a = - (8) 

where e is the dielectric constant of silicon dioxide, and f Q is the permittivity of 
free space. The initial linear density, N 0 , depends on particle energy and can be 
found from stopping power tables (Ref 10:310). The value of the Gaussian radius, 
b, was taken from Reference 7 where it was used as a variable meter to fit 
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this recombination model to experimental data. The value was found to be 3.5 ' 

nanometers (Ref 7:2697). J 

It should be pointed out that this model is only valid for situations where the | 

electron-hole pair separation distance is short enough that the pairs lose their [ 

identity. In this case electrons readily recombine with holes from different pairs. 

This is generally true of the particle tracks left by heavy charged particles, 

because of their high linear energy transfer. For electrons and gamma-rays, the ; 

pair separation distance is large enough that recombination is pretty much res¬ 
tricted to electron-holes of the same pair. The model developed to handle these 1 

cases is called the geminate model (Ref 1:27). j 

In general, the yield fraction f y depends on three factors. First, it depends on 
the energy of the particle. This is because the initial linear density, N 0 in Eq (3), ! 

generally decreases with increasing particle energy, at least in the ranges of \ 

interest here. The lower the initial density, the less probability there is that 
recombination will occur, thus the yield fraction increases. Second, yield fraction 
also depends on the strength of the electric field. The stronger the field, the more 
quickly the electron and hole columns separate, decreasing the time available for 
recombination. Thus, yield fraction increases with increasing electric field. Third, 
the yield fraction depends on the angle of incidence. This can be understood by 
considering the two extremes. When the particle track is perpendicular to the 
field, the two columns are separated rapidly, which increases yield fraction. When 
the particle track is parallel to the field, the columns are pulled through each 
other, enhancing recombination, and reducing the yield fraction. Thus, yield frac¬ 
tion increases with increasing angle of incidence. 

This increase in yield fraction with particle energy, electric field and angle of 
incidence is the cause of the general characteristics of the curves seen in Figures 3 
through 7. In these figures, the change in threshold voltage also increases with 
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particle energy, electric field, and (except for particle energies below 6 MeV) angle 
of incidence. As will be seen below, the change in threshold voltage is directly 
proportional to the charge accumulated at the interface of the oxide and thus the 


Bulk Phenomena 

After escaping initial recombination the remaining electrons are rapidly swept 
from the oxide by the applied gate field. This process takes a few picoseconds to 
complete. Two assumptions are made about this stage. First, no electrons are 
trapped within the oxide. Second, having escaped columnar recombination, the 
electrons do not recombine further with holes previously trapped in the oxide. 

After the electrons are swept from the device, the holes transport toward the 
interface much more gradually. Assume all the holes move with the same velo 
city, v = fj, + E. For an oxide of length 280 nanometers and an electric field of 
1 X 10 3 Y/cm at room temperature, it takes a hole 1.40 X 10 -5 seconds to cross 
the oxide. This time is short enough to justify ignoring any transient effects 
caused by transport. To simplify the analysis then, all the holes produced by an 
incremental dose, AD (that escape recombination), are assumed to instantaneously 
transport to the interface. Bulk trapping within the oxide is ignored. 


Interface Phenomena 

.As the holes pass through the interface region, a certain fraction, f T , are 
caught in deep hole traps within about 10 nanometers of the interface. The frac¬ 
tion trapped depends on the temperature and electric field, but is mostly a result 
of the processing techniques used during manufacture. The interface region has 
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high local strain and a deficiency of oxygen atoms. This gives rise to a large 
number of Si-Si bonds instead of the usual Si-O-Si bonds found in silicon dioxide. 
When a hole passes one of these bonds, it can break the bond and recombine with 
one of the bonding electrons. This results in a Si atom retaining the remaining 
electron from the broken bond, and the positive charge residing with the other Si 
atom (Ref 1:34). 

Since the trapping fraction depends primarily on the manufacturing process, 
it can vary quite radically between different MOSFETs. Generally, f T can be as 
little as 1-2 percent for hard oxides, 10-20 percent for good quality commercial 
oxides, and as much as 50-70 percent in very soft commercial oxides (Ref 1:33). 

The spatial dependence of the charge trapped at the interface is also depen¬ 
dent on the relative hardness of the device. Hard oxides have distributions that 
fall off exponentially as one moves away from the interface into the oxide. On the 
other hand, the softer oxides have distributions that remain relatively constant 
(Ref 11:1207). 

These trapped holes can persist from a few hours to years. Annealing does 
take place, in which electrons from the Si semiconductor tunnel into the oxide and 
recombine with the trapped holes. However, since the time frame for annealing is 
considerably longer than the experiment, it is ignored in this model. 

This model assumes that a certain fraction, f-p, of the holes escaping recombi¬ 
nation will be trapped in a constant charge distribution, p, extending a distance X 
from the interface. The incremental charge density due to an incremental dose is 
then 


-V ■ 


KgADIf y f T 


(9) 


in the region near the interface and 0 elsewhere. X is presumably about 10 nm 
(Ref 1:33). The trapping fraction is uncertain, but since these devices are reported 




16 








as non-radiation hardened (Ref 3:4393), it is probably greater than 10 percent. 
This quantity is left undetermined in the analysis, but is assumed constant. 

Radiation-induced interface traps exist very close to the interface, and within 
the Si band gap. Their occupancy is determined by the position of the Fermi 
level at the interface, and can be negative (usually with n-channel devices), posi¬ 
tive (usually with p-channel devices), or neutral. In either case, the presence of 
these trapped charges results in shifts to the threshold voltage. At the present 
time, the precise mechanisms responsible for the formation of these states are not 
well understood, and a subject of considerable debate (Ref 1:37). These states are 
difficult to model. Since the deep hole traps are dominant (especially at the lower 
dose levels applicable here), interface traps are ignored in this analysis. 


Threshold! Voltage vs Dose 

In this section, an analytical expression relating the change in threshold vol¬ 
tage to the total ionizing dose is found, so that the curves in Figures 5 through 7 
can be simulated. This expression will include the dependence yield fraction has 
on electric field and incident angle. 

The change in threshold voltage for arbitrary charge density is given by (Ref 
12:199) 

AV, h -/^W<ix (10) 

0 t 

where 

1 = the length of the oxide (cm) 
p(x) = density of holes (cm -3 ) 
e = permittivity of Si0 2 


.Assuming the charge density is constant between 1 — X and 1, and zero elsewhere. 












the threshold voltage change becomes 


AV lb = i2^SI 


dAV th -^- 2XI I ( 

th 2e 

where dAV th is the incremental change in threshold voltage resulting from an 
incremental change in charge density. Substituting the value for A p given by I 
(9) into Eq (10) gives 

dAV lh - K & f y f T l (X -~ 2I) AD ( 

In Eq (13), everything is assumed constant except the yield fraction, f v , which 
depends on the electric field, which in turn depends on the charge distribution. 

We next derive an expression relating the internal electric field to p and th 
AV th . Figure 10 shows a schematic cross section of a MOSFET, with an areal 
density Q x - located at arbitrary x' within the oxide. This charge induces a cha 
density Q, at the gate, and Q ; at the semiconductor. Ignoring any difference in 
work function between the metal gate and semiconductor, it can be shown that 


Using Gauss’ Law, the internal electric field due to charge density Q x . becomes 
E(x) = I 


for 0 < x < x' and 




Figure 10. Schematic Cross Section of MOSFET With Arbitrary Areal Chs 
Density Q x . at x' 

E(x) - ^-Q„. I 

for x' < x < 1. 

Eqs (16) and (17) can be generalized to an arbitrary charge distribution p(: 
as shown in Figure 11. Taking into account the field due to the external gate \ 
tage and letting Q x = p(x' )dx', the electric field at x becomes 

E(x) = /^L^lD-Hx' - Mx') , ix , + i , 

0 f x f ‘ 

Expanding the second term on the right and rearranging gives 


















dAV, 



Integrating Eq (22), and rearranging, results in an equation relating the threshold 
voltage change and the dose: 




£ 
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HI. Results and Analysis 

The previous section developed the model used to simulate the effect charged 
particle radiation has on MOSFET devices. This section describes the results of 
applying this model to the actual experiment described in the introduction and 
background section. 

Figure 12 shows the yield fraction as a function of particle energy for angles 
of incidence of 0, 45, and 80 degrees using the columnar recombination model. An 
electric field of 1.79 X 10 s V/cm was used. This field was found by dividing the 
gate voltage (5 Volts) by the thickness of the oxide (280 nanometers). The fraction 
yield increases with both particle energy and angle of incidence, and the charac¬ 
teristics of these curves are similar to those shown in Figures 3 and 4 except for 
the data anomaly. 

In finding a functional relationship between the change in threshold voltage 
and the dose (accounting for oxide field modulation due to hole trapping near at 
the interface) it was assumed that the fraction yield was a linear function of the 
electric field within the oxide. Figure 13 shows the fraction yield as a function of 
electric field for 4 MeV protons with angles of incidence of 45 and 80 degrees. 

This figure was found by using the columnar recombination model. .As can be 
seen, the fraction yield increases with both electric field and angle of incidence, 
but the curves are not linear throughout the range of interest. To properly use 
Eq (23). a linear approximation was taken from 0.1 NIY/cm to 0.25 MV/cm (the 
general range of the oxide electric field during the irradiation). The coefficients in 
Eq (21) are found to be a* = 4.11 X 10 - ' cm/V and bj = .0539 for an angle of 45 
degrees, and a, = 4.88 X 10 -7 cm/V and b, = .0650 for an angle of 80 degrees. 
hVing these values in Eq (23) with X = 10 nm, 1 = 280 nm. and \ g = 5 Volts, 
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Figure 12. Fraction Yield as a Function of Particle Energy and Angle of 
cidence. 

AV th was found as a function of dose for trapping fractions of .20, and .10. Th 
results are shown in Figure 13. As expected, the threshold voltage rises with do 
angle of incidence, and trapping fraction. However, the sought after saturation 
effect is missing. 

The reason the saturation effect is absent is straightforward. If all of the 
charge is assumed to accumulate within 10 nm of the interface, the charge givin 
rise to a threshold voltage change equivalent to the gate voltage will (from Eq 
(20)) only decrease the internal field to 0.982 of its original value. This is gen¬ 
erally true for charges at the interface—they modulate the internal electric field 
very little, certainly not enough to cause th*' ^turation effect seen in Figures ■’> 
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these parameters, it is likely that a significant proportion of charge is trapped in 
the bulk of the oxide (Ref 13:3940). 

The implication of large bulk trapping can be seen from Figure 15. This 
figure shows the electric field as a function of distance for a charge density distri¬ 
buted evenly throughout the oxide. The charge density is assumed large enough 
to give rise to a threshold voltage change equal to the applied gate voltage. From 
Eqs (11) and (19) the internal electric field is given as 

E(x) = 2V g -l (24) 

so E(0) = 0 and E(l) = 2V g /l. The distributed charge causes greater modulation 
throughout the bulk of the oxide, than that charge accumulated at the surface. 

To the left of the midpoint in Figure 15 the electric field falls significantly, 
causing a reduction in yield fraction. On the other hand, the field rises 
significantly to the right of the midpoint, which increases the yield fraction. 

From Figure 13. the yield fraction seems to slope off with decreasing electric field, 
so the decline in fraction yield in the left hand region may dominate the increase 
in the right hand side, resulting in an overall reduction in fraction yield. Whether 
this reduction is great enough to cause the saturation effect remains questionable. 

It could well be a contributor, but there are other things to consider. 

For instance, it was assumed earlier that the electrons, having escaped colum¬ 
nar recombination, were swept from the oxide without further recombination. 

F'or low hole densities and high electric fields this assumption is certainly valid. 
However, the cross section for electron-hole recombination is proportional to E -i 2 
(Ref 14:3203). The reduced field in Figure 15 could result in significant recombi¬ 
nation in the region to the left of midpoint. This could be another source of the 
saturation effect, especially if the field reaches zero within the oxide. Then the 
electrons would stop and recombine with a hole (Ref 16:1522-1523). 
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Figure 15. Electric Field Within the Oxide from an Evenly Distributed Charge 
Density 

A third assumption made in this model was that the fraction of holes 
trapped at the interface was constant. However, it is known that the hole trap¬ 
ping cross section is dependent on the electric field and varies with E -1 2 (Ref 
15:1192). The increasing field at the interface, as seen in Figure 15 may result in 
a decline in trapping fraction, that could also be a contributor to the saturation 
effect. 

These three field modulated mechanisms—yield fraction, electron recombina¬ 
tion, and hole trapping may all contribute to the saturation effects seen in Figures 
5-7. To verify this will require a more refined model that not only takes these 
mechanisms into account, but also accounts for the wide iispersion inherent in the 
hole transport. This will be required to characterize the hole buildup in the bulk 








W. Conclusion 


In conclusion, this simple model has shown that the electric field reduction 
due to holes trapped at the interface cannot cause the saturation effect. Charges 
trapped close to the interface cannot reduce the internal electric field to any great 
degree--certainly not to the extent required to effect the yield fraction of holes 
escaping columnar recombination. However, similar p-channel and n-channel dam¬ 
age sensitivities indicate that considerable charge is trapped within the bulk of the 
oxide. Charge distributed within the oxide bulk modulates the electric field con¬ 
siderably. This modulated electric field could contribute to the saturation effect in 
three ways— reduction of the yield fraction throughout the oxide, increase in elec¬ 
tron recombination, reduction of hole trapping at the interface. 

To verify this, the present model would have to be refined in a number of 
ways. First, the dispersive nature of hole transport would have to be considered 
so that reasonable estimates of the trapped hole distribution (both at the interface 
and within the bulk of the oxide) can be made. Second. Eqs (10) and (19) would 
have to be modified to numerically find the threshold voltage and electric field for 
an arbitrary charge density. Third, the nonlinearity between yield fraction and 
electric field would have to be accounted for, so that overall yield fraction due to 
a widely varying internal electric field may be found. Fourth, the electron tran¬ 
sport and electric field dependent recombination cross section would have to be 
introduced. Finally, the hole capture cross section dependence on electric field 
would have to be taken into account. Introducing these refinements into the 
model should allow determination of the primary cause of the saturation effect 
and explain the data anomaly. 
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Appendix: Columnar Recombination Computer Code 


This appendix contains the code used to find the yield fraction of electro 

and holes escaping initail recombination. The following algorithm is used: 

Step 1. Input the initial linear density, the mobilities, the diffusion coefficie 
the recombination coefficient, the electric field, the angle of inci 
incidence, the column radius, and the length of the oxide (lines 22-: 

Step 2. Set up the x and y grid. In this case the maximum values in the x 
y directions and the number of grid spaces are read in. The progra 
set up to take advantage or the svmmetrv in the v direction (lines : 
45). 

Step 3. Input and initialize time related elements. These include the time i 
ment and the ratio of time increments between each step. This rat 
inserted to help shorten the run times. The time variable is initiali 
and the maximum time calculated (lines 51-53). 

Step 4. Set up boundary conditions on the grid. In this routine the bound; 

are assumed large enough that the electrons and holes separate well 
before the columns reach the boundary time so the electron and ho 
densities at xmax. xmin. and ymax are all set to zero (lines 58-73). 

Step 5. Determine the linear density of particle at time equal to zero (lines 
124). 

Step 6. Determine the yield fraction for each time step (lines 132-216). 

Step 7. Determine new time and time ratio (lines 133-134). 

Step 8. Set new boundary condition along ymin (lines 139-142) 


Step 9. Determine the new density at each point on grid (lines 144-183 


Step 10. 
Step 11. 
Step 12. 
Step 13. 


Kind the diffusion term in Eq (2) (lines 151-160). 
Kind the drift term in Hq (2) (lines 164-170). 


Kind the recombination term in Eq (2) (line 171). 

bind the time rate of change of the holes electrons (lines 1 








Step 15. Find the number of charge carriers per unit length and thus th 

new yield (lines 200-212). 

Step 16. Print out the time,electron and hole density, and the yield (lii 

213). 

Appendix of variables. 

alpha - recombination coefficient in cm3/sec 
b - Gaussian radius in nanometers 
bcm - Gaussian radius in centimeters 

check - variable used to determine whether the point is outside the column or 
delx - the length of each element in the grid in the x- direction in nanometers 
dely - the length of each element in the grid in the y- direction in nanometers 
delt - the time increment in seconds 

derele(i.j) - time rate of change of electron charge density in /cm3*s 

derhol(i.j) - time rate of change of hole charge density in /cm3*s 

diffn - diffusion term in the electron density rate of change equation 

diffp - diffusion term in the hole density rate of change equation 

dndx - first derivative with respect to x of electron density 

dpdx - first derivative with respect to x of hole den density 

dpi us - diffusion coefficient of holes in cm2/s 

driftn - drift term in the electron density rate of change equation 

driftp - drift term in the hole density rate of change equation 

dneg - rhp diffusion coefficient of the electrons in cm2/s 

d2ndx2 - second derivative with respect to x of the electron density 

d2ndy2 - second derivative with respect to y of the electron density 

d2pdx2 - second derivative with respect to x of the electron density 

d2pdy2 - second derivative with respect to v of the hole density 






e - the electric fiefield applied across the oxide in Volts/cm 
eleden(i,j) - the electron charge density in charges /cm3 

eledenx - the density of electrons along the x axis 

eledeny - the density of electrons along the y axis 

elenum - the number of electrons along the y axis in charge/cm 

elenuml - the number of electrons above the y axis in charge/cm 
elenumo - the total number of electrons equal to elenum+elenuml*2 
hx - thethe length of each element in cm in the x direction 
hy - the length of each element in cm in the y direction 
holden(i,j) - the hole charge density in charges /cm3 
holdenx - the hole density along the x axis in charges/cm3 
holdeny - holthe hole density along the y axis in charges/cm3 
holnum - the number of holes along the y axis in charges/cm 
holnuml - number of holes above the y axis in charges/cm 
holnumo - total number of holes equal to holnum+holnuml*2 
length - length of oxide in cm 

m - number of elements between upper and lower boundaries in the grid 
direction 

midx - midpoint of the grid in the x direction 
midv - midpoint of the grid in the y direction 
muneg - mobility of electrons in cm2/V*sec 


n - number of elements between the upper and lower boundaries of the i 
y direction 

no - initial line density in charges/cm 

recom - recombination term in the density rate of change equation 

thetad - angle of incidence between the particle track and the electric fie 
degrees 

thetar - angle of incidence in radians 
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This code is designed to find the concentration of holes and elec- 
trons(in charges/cm**3) and the yield fraction of holes escaping 
recombination aue to interaction of charged particles (alphas and 
protons) with silicon dioxide. 

implicit none 

real pi,no,muneg,muplus,dneg,dplus,alpha,e,thetad,thetar 
real b,bcm,xmax,xmin,ymax,ymin,delx,dely,x(201),y(102) 
real delt,eleden(201,102),holden(201,102),check,yield 
real derele(201,102),derhol(201,102),hx,hy,tmax,time,length 
real elenum.holnum 

real elenumo,holnumo,eledenx(201,100),holdenx(201.100) 

real eledeny(l01,100),holdeny(101.100),elenuml,holnuml,d2ndx2 

real d2ndy2,diffn,dndx,driftn,recom,tratio 

real d2pdx2.d2pdy2,diffp,driftp.dpdx 

integer n,m,midx,midy,i,j 

open(unit=l,file=’rec.dat\status=’old’.READONLY) 


This section initializes a number of the constants. 

read(l,*)no,muneg,muplus,dneg,dplus,alpha 

read(l,*)e,thetad,b,length 

pi=3.1416 

thetar=thetad*pi/180 


: This section sets up the x and y grid. 

read(l,*)xmax,xmin,ymax,ymin.m.n 

bcm=b*lE-07 

midx=m/2+l 

midv=2 

delx=f xmax-xminl/m 
dely=(ymax-ymin)/n 
hx=delx*lE-07 
hy=dely*lE-07 
x(li=xmin 
y(l)=ymin 
do 10 i=2.m+l 
x(i)=x(i-l)+delx 
10 continue 

do 20 j=2,n+l 
y(j)=y(j-l)+dely 


This section inputs the time related elements such f 
increment, and the maximum time considered. 

read(l,*)delt.tratio 

tmax=(length/cos(thetar))/(muneg*cos(thetar)*e) 

time=0.0 
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0055 

0056 

c 

This section inputs the boundary conditions for the charge density 

0057 

c 

at xmin and xmax. 

0058 

0059 


do 30 j=2,n+l 

0060 


eledenl l,j)=0.0 

0061 


eleden(n+l,j)=0.0 

0062 


holden(l,j)=0.0 

0063 


holden(n-f l,j)=0.0 

0064 

30 

continue 

0065 

0066 

0067 

c 

This section inputs the boundary conditions for the charge density 

0068 

c 

at ymax. 

0069 

0070 


do 40 i=2,m 

0071 


eleden(i,n+l)=0.0 

0072 


holden(l,n+l)=0.0 

0073 

40 

continue 

0074 

0075 

0076 

c 

This section initializes the charge densities at time equal zero. 

0077 

0078 


do 50 i=2,m 

0079 


do 60 j=2,n 

0080 


check=-(x(i)**2+y(j)**2)/b**2 

0081 


if(check.lt.-100.00)tnen 

0082 


eleden(i,j)=0.0 

0083 


holden(i,j)=0.0 

0084 


else 

0085 


eleden(i,j)=(no/(pi*bcm**2))*exp(-(x(i)**2+y(j)**2)/b**2) 

0086 


holden(i,j)=eleden(i,j) 

0087 


endif 

0088 

60 

continue 

0089 

50 

continue 

0090 

0091 

0092 

c 

This section finds the charge density along the x and v axis at 

0093 

c 

time equal to zero. 

0094 

0095 


do 70 i=l,m + l 

0096 


eledenx(i.l)=eleden(i,midy) 

0097 


holdenx(i,l)=holden(i,midy) 

0098 

70 

continue 

0099 


do 80 j=l,n+l 

0100 


eledeny(j,l)=eleden(midx.j) 

0101 


holdeny(j,l)=holden(midx.j) 

0102 

80 

continue 

0103 

0104 

0105 

c 

This section determines the total number of charge carriers(per 

0106 

c 

unit length) at time equal to zero. 

0107 

0108 


elenum=0.0 
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0109 holnum=0.0 

0110 elenuml=0.0 

0111 holnuml=0.0 

0112 do 90 i=2,m 

0113 do 300 j=3,n 

0114 elenuml=elenuml+eleden(i,j)*(hx*hy) 

0115 holnuml=holnuml+holden(i,j)*(hx*hy) 

0116 100 continue 

0117 90 continue 

0118 do 110 i=2,m 

0119 elenum=elenum+eleden(i,2)*hx*hy 

0120 holnum=holnum+holden(i,2)*hx*hy 

0121 110 continue 

0122 elenumo=elenum+elenuml*2 

0123 holnumo=holnum+holnuml*2 

0124 yield=holnumo*delt 

0125 print *,time,elenumo,holnumo,yield 

0126 print*,'max time=’.tmax 

0127 

0128 

0129 c This section determines the time derivative of the charge density 

0130 c the new charge density and the total charge as a function of time. 

0131 

0132 120 if(time.gt.tmax)go to 130 

0133 time=time+delt 

0134 delt=delt*tratio 

0135 elenum=0.0 

0136 holnum=0.0 

0137 elenuml=0.0 

0138 holnuml=0.0 

0139 do 140 i = l,m + l 

0140 eleden(i,l)=eleden(i,3) 

0141 holden(i,l)=eleden(i,3) 

0142 140 continue 

0143 

0144 do 150 i=2,m 

0145 do 160 j=2,n 

0146 
0147 

0148 c This section finds the diffusion term for the electrons from finite 
0149 c difference form. 

0150 

0151 d2ndx2=(eleden(i+l,j)-2*eleden(i,j)+eleden(i-l ,j))/hx* *2 

0152 d2ndy2=(eleden(i,j + l)-2*eleden(i.j)+eleden(i.j-l ))/hy**2 

0153 diffn=dneg*(d2ndx2+d2ndv2) 

0154 

0155 c This section finds the diffusion term for the holes from finite 
0156 c difference form. 

0157 

0158 d2pdx2=(hnlden(i + l ,j)-2*holden(i.j)+holden(i-1 ,j))/hx**2 

0159 d2pdv2=(holden(i.j + l )-2*holden(i.j )-*-holden(i.j- l))/hv**2 

0160 diffp=dplus*(d2pdx2+d2pdv2) 

0161 

0162 <• This section finds the drift term for the electr 








dndx=(eleden(i+l ,j )-eleden(i-l ,j ))/(hx*2) 
driftn=muneg*e*sin(thetar)*dndx 

This section finds the drift term for the holes. 

dpdx=(holden(i+l,j)-holden(i-l,j))/(hx*2) 

driftp=muplus*e*sin(thetar)*dpdx 

This section finds the recombination term for both holes and electrons. 

recom=-alpha*eleden(i,j)*holden(i,j) 

This section then finds the time rate of change of the electron and 
hole density for each x and y. 

derele(i,j)=diffn+driftn+recom 

derhol(i,j)=diffp+driftp+recom 

continue 

continue 

This section finds the new charge density at each point for 
new time. 

do 170 i=2.m 
do 180 j=2.n 

eleden(i,j)=eleden(i.j)-Hlerele(i.j)*de)t 

holden(i.j)=holden(i.j)+derhol(i,j)*delt 

if(eleden(i,j).lt.0.0)then 

eleden(i,j)=0.0 

endif 

continue 


This section finds the number of charge carriers(per unit length). 

do 190 i=2.m 
do 200 j=3.n 

elenuml =elenuml+eleden(i.j)*(hx*hy) 
holnuml = holnuml+holden(i.j)*(hx*hy) 
continue 
continue 
do 210 i=2.m 

elenum=elenum+eleden(i,2)*hx*hy 

holnum=holnum+holden(i.2)*hx*hy 

continue 

elenum=elenuml*2+elenum 
holnum=holnuml*2+holnum 
yield=yield+holnum*delt 
print*,time.elenum.holnum, yield 
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19. Abstract 

A simple model to characterize the effects of heavy charged particles 
on metal-oxide semiconductor field-effect transistors (MOSFET) is presented. 
The model is applied to experimental results provided by the Air Force Wea¬ 
pons Laboratory, and an attempt made to simulate saturation phenomena ob¬ 
served in the threshold voltage change. The model assumes all trapped holes 
are within a few nanometers of the oxide-semiconductor interface, and takes 
into account the resultant counter electric field, and its effect on the frac¬ 
tion escaping recombination. An equation relating threshold voltage change as 
a function of dose is derived and used, but does not duplicate the saturation 
effect. This is because the charge trapped near the interface reduces the in¬ 
ternal oxide electric field very little. However, similarities in p-channel 
and n-channel MOSFET damage sensitivities indicate considerable bulk charge 
trapping. This would modulate the internal electric field considerably. To 
determine whether this accounts for saturation will necessitate refinement 
of the model, taking hole transport and bulk trapping into account. 
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